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Abstract 

We propose a microscopic traffic model where the update velocity is determined by the decel¬ 
eration capacity and response time. It is found that there is a class of collisions that cannot be 
distinguished by simply comparing the stop positions. The model generates the safe, comfortable, 
and efficient traffic flow in numerical simulations with the reasonable values of the parameters, and 
this is analytically supported. Our approach provides a new perspective in modeling the traffic-flow 
safety and the perturbing situations like lane change. 
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I. INTRODUCTION 


Modeling of traffic flow lias been an intensive research topic for more than a half century 
in the engineering and science communities [l- 


reviews 


12 


111 ], of which results are summarized in the 


15]. In the progress of information technology, the various traffic models are 


required for the control and/or the automation of traffic flow. The needs is basically the 
credible modeling of safety and mobility, the two categorical but conflicting goals in driving. 
Thus the natural driving behaviors have been modeled, for example, the more (less) accel¬ 
eration for the larger (smaller) spacing. However, this is usually based on the trial-function 


approach, as criticized in Ref 


12 ]- 


There were a few seminal works that do not use trial function; one is Gipps (collision- 
avoidance) model [7, 12] and another one is Nagel-Schreckenberg (minimal collision-free) 


model 


13]. In Nagel-Schreckenberg model, any magnitude of deceleration is applied 


when required to prevent a collision. This means the deceleration capacity is actually un¬ 
bounded (the so-called intelligent-braking-behavior suggested in [ll] also belongs to this 
case). Meanwhile, in Gipps model, a collision-avoidance in bounded deceleration capacity 
was suggested. We consider this approach is more physical, and thus we adopt it in the 
present work. 

In this work, we examine how the collision can be understood in the physical constraints 
of the deceleration capacity and the response time. We propose a microscopic traffic model 
where the update velocity is determined in the safety criterion by the constraints. It is 
found that there is a class of collisions that cannot be identified by the usual safety crite¬ 
rion comparing the emergency-stop positions. The resultant model generates in numerical 
test the practically appealing traffic flow of safety, efficiency, and comfort with the reason¬ 
able parameter values, and this is analytically supported. Our model also provides a new 
perspective in modeling traffic-flow safety and the perturbing situations like lane change. 


II. MODELING 

Since the safety and mobility are conflicting to each other, a compromise between them 
is necessary. A natural one is the condition where driver can marginally avoid collision 
against the leader’s full stop. Let x n (t ) and v n (t) be the (front-end) position and velocity at 
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FIG. 1. (Color online) Two types of safety criterion in emergency when vehicles decelerate by 
their own braking capacities till stop. Follower’s deceleration as the response to the emergency is 
delayed by response time r n (the dotted segment from t and t + r n is not the part of emergency). 
The curvature of trajectory is given by the associated braking capacity. In criterion (a), the two 
stop positions are compared to tell a collision. In criterion (b), the blue and black solid trajectories 
meet though the follower’s stop position does not exceed the leader’s (this is possible only when 
the follower’s braking capability is stronger than the leader’s). Note this kind of collisions cannot 
be distinguished by comparing the stop positions (see the red dashed curve). 

time t, respectively, of vehicle n. The safety criterion is asking, in the presence of response 
time r n , what is the marginal x n {t + r n ) and v n {t + r n ) that does not result in a collision 
with the maximum deceleration D n from t + r n , if the front vehicle at x n+ \ (t) and v n+ \ (t) 
begins to decelerate with its maximum deceleration D n+ 1 from t to stop. In short, this asks 
whether the worst situation is manageable in the physical constraint of braking capability 
and response time. Obviously, such a worst case may not happen. But it is necessary to 
check whether the follower can keep safe in that situation with its braking capacity and 
response time. 

One of the safety criteria is shown in Fig. [0(a), where x n (t+r n ) and v n (t+r n ) are adjusted 
so that the two trajectories become tangential as the two vehicles stop. It basically compares 
the two vehicles’ stop positions to tell a collision. This is same to that considered in Gipps 
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model [7] believed so far to provide a safe enough dynamics, ffere, we point out that this 
safety criterion only is incomplete. This is because there is the other kind of collisions that 
cannot be discerned by comparing the stop positions, as follows. 

The other kind is shown in FigQJb), which is possible only when D n > D n+ \. In this 
case, since the follower’s trajectory is bent stronger than the leader’s, the match of the stop 
positions (see the red-dashed curve) necessarily brings about a collision before stop. This 
collision is, however, not distinguished by simply comparing the stop positions. It is thus 
necessary to reconsider the configuration at t + r n . The follower’s blue solid trajectory in 
Fig. [IJb) is the alternative, which is adjusted to be tangential to the trajectory of the leader 
still in move. We remark that, even for D n > D n+l , there is a situation where the criterion 
of Fig. [0(a) should still apply, for example, if the follower is not so close to the leader at 
time t. 

x n {t+T n ) and v n {t + r n ) are related by position-update rule. When the scheme of constant 
acceleration between responses is used, the position update reads 

X n {t + T n ) = X n (t) + y [v n {t) + V n (t + T n )\ . (1) 

Considering this in the two tangential conditions explained above, as the two marginal 
velocities at t + r n , one can obtain 

Vn(t + r n ) = -yy + ^ + D n S n (t) - T n V n (t ) + 

^ - AD n (2 s n (t) + r n Av n (t)) , 

where A D n = D n+1 - D n: A v n (t) = v n+1 (t) - v n (t), and s n (t) = x n+ 1 (t) - x n (t) - L n+1 for 
the leading vehicle’s length L n+ 

V ^(t + T n ) is enough to tell a collision when D n < D n+ 1, while it is not when D n > D n+1 . 
Thus in the latter case, one of v s n (t + r n ) or v^(t + r n ) should be selected depending on 
situations. Although Fig. [0(b) shows a situation where v^(t + r n ) should be selected, this is 
not always the case. If the follower is not so close to the leader, one can easily argue that 
v n(t + T n ) is instead the proper choice. This way, considering a few conditions, one knows 


r n 


v n(t + T n ) — V n+ 1 (t) — (D n + D n+ 1 ) + 


'r n A D 


J n +1 


H) 


D 


n +1 


( 2 ) 
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the candidate of the update velocity is given by 


,cand 


(t+T n ) = < 


v n(t + T n) if D n > D n+ 1 , v d {t + T n ) is real, and r n + < 15 ^, 

v n(t + T n ) else if v„(t + T n ) is real, 

v n if) — T n D+ otherwise, 


( 3 ) 

where the last case is introduced to cover such a situation allowing no physically meaningful 
v^(t + r n ) and v s n {t + r n ). For example, a careless cutting-in may bring it about. This means 
there is no way to avoid a collision if the cutter brakes maximally till stops. D+ is introduced 
as an indicator of such an emergency that requires a deceleration larger than D n , which is 
not possible by the definition of D n . In order to simply cover that situation in the model, 
one may assign a value larger than D n to D+. We anticipate that the third case is crucial 
in modeling (un)tolerable perturbations. 

v™ nd (t+T n ) in Eq. (J3j) is still the candidate velocity of the next step because its realizability 
from the current velocity is not taken into account yet. The realizability is determined in 
the vehicular performance represented by the deceleration and acceleration capacities. Thus 
when an acceleration capacity A n is additionally introduced, the realizability corresponds 
to D n < (t£ and (t + r n ) —v n (t))/r n < An”, named as “mechanical restriction” jlfij] . When 
a velocity change exceeding this range is required, only —r n D n or r n A n is possible by the 
definition of D n and A n . Finally, considering the traffic regulation also, we arrive at 


v n (t + Tn) = max{w n (f) - T n D n , 0, min{v n (f) + T n A n , u max , < and (f + r n )}}, (4) 


where zero stands for the directionality and n max is a speed limit. This gives the velocity 
update with Eqs. (J 2 ]) and d3J), and then position update is followed in Eq. (jT]). The update 
rule is applied in parallel to all vehicles in system. 


III. MANAGEABILITY AND POTENTIAL COLLISION 

Before analyzing our new model, we discuss a few implications of Eq. (4). The interest 
is in the case when 

< and (^ + Tn) > V n (t + T n ). (5) 

We below call a vehicle holding Eq. ([5j) manageable at time t. The inequality says the 
realization of v™ nd (t + r n ) is possible in the braking capacity D n . This again implies, if 
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a deceleration is required for safety, it is realizable. Interestingly, the manageability at t 
lasts thereafter unless perturbed later. This is attributed to the fact that v^ d (t + r n ) are 
constructed in a way to keep safety against the leader’s worst behavior; the consecutive 
maximal braking to stop (if already stopped, it is assumed not to move). Thus once all the 
vehicles in a system are manageable, it lasts forever and the traffic flow remains collision- 
free, as long as no perturbation is applied. A closed system composed of vehicles initially 
at rest is such an example. 

As a perturbation, one may consider the insertion of a vehicle into a gap between two 
vehicles. When the insertion takes place, it is reasonable to examine the manageability of 
the follower and that of the new comer. One may call it manageable insertion when the two 
vehicles are manageable at that instant. We consider the notion of manageable insertion is 
crucial in modeling on-ramp and/or lane-change. Also, a strategy to the dilemma zone by 
the traffic signal can be examined as considering the insertion of a standing object. This 
way, the manageability may shed a light in designing and/or modeling the (un)tolerable 
traffic perturbations. 

From the other perspective, the non-manageability [violation of Eq. (jSJ)] can give a mea¬ 
sure for safety indicating a possible collision. The situation of non-manageability results in 
collision if the leader really applies the maximal brake to stop. Thus the statistics on the 
non-manageability can be a reasonable measure for the traffic-flow safety. Note the flow 
including non-manageable configurations does not necessarily result in collision. Therefore, 
the flow without a collision can be regarded as dangerous in our approach. We believe this 
viewpoint should be applied to the real traffic. It is necessary to discern a traffic flow with 
potential collisions so as to prevent a traffic accident in advance. 


IV. NUMERICAL RESULT AND ANALYTIC SUPPORT 


The present model is a consequence of the physical meaning of D n and r n . We thus 
examine the flow property while varying them in numerical study. For the other model 
parameters, we use the typical values; u max = 110 km/h, A n = 1.5 m/s 2 , and L n = 7.5 
m [l6, 17]. A system of randomly distributed vehicles initially at rest in the 100-km-long 
circular road is tested for various vehicular density p (the total number of vehicles divided 
by the road length). For D n , a random value out of the interval (rrio — + wjj) is 
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FIG. 2. (Color online) Three phases, (a) The location-speed snapshots for r = 1 sec. of hS flow 
(p = 30 veh/km, mp = 5 m/s 2 ), fS flow (p = 68, mo = 5), and jam phase (p = 97, mn = 10). (b) 
The phase diagram in the parameter space (p,mu) for r = 1. Dotted (solid) lines are numerically 
(analytically) obtained phase boundaries. 


assigned for various and fixed % = 1.5 m/s 2 (wp -value turns out not to change the 
results qualitatively). Below, we will use r n = r for all n for a simplicity of the numerical 
implementation. 

We observe that there emerge three kinds of steady state traffic flow depending on mo, t, 
and p. Figure 0(a) shows the three flows with the position-velocity snapshots, where each dot 
represents a car. One of the snapshots (circles) shows an almost flat velocity profile, and we 
name it “homogeneous steady” (hS) flow. Another snapshot (squares) exhibits fluctuating 
velocities, named as “fluctuating steady” (fS) flow. Meanwhile, the last one (crosses) shows 
a traffic jam (J) where vehicles can hardly move. Figure E](b) shows the phase diagram on 
p-rrto plane for r = 1 sec., where the dotted (solid) curves are the numerically (analytically) 
obtained phase boundaries. We observe that the boundary between hS and fS is robust, 
while the region for J phase depends on the initial configuration. In the following, we 
analytically argue that the observation above is the intrinsic feature of our model. 

The hS flow is a homogeneous-velocity solution (HVS) of the model. Considering a 
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constant velocity v in Eq. (JHJ) for all n, one can find (see Appendix A) 


P 


-i 


V 

TT] 


T 2 T] 


/(^) 


for — < 1 , 

TTj — 7 

for - 2 - > 1 

TT] 


( 6 ) 


for rj 1 = (m D -w D ) 1 -(m D + w D ) 1 , where f(z) = (z - z 2 y/2) + f* /z l/2y p(y) dy 
for p(y) the probability distribution of the random variable y standing for r)(l/D n+ i — 1 /D n ) 
[Eq. (IA22 h is the details of f(z) for the D n s we use]. As a general property indifferent to the 
statistics of D n s, f(z) is increasing and convex downward, and f(z) —> z for z —> 1 while 
f(z) ~ z 2 for z 1 [see the two limiting behaviors of the red curve (and also data points) 
in the upper-right part of Fig. |3](a)]. 

The solid curve in Fig. [3](a) is Eq. (1A22I) [a realization of Eq. (EJ) for the D n s we use]. 
The numerical data for hS are perfectly on it above v/rr) = 1. Therein, the circles and 
bars are, respectively, the velocity averages and fluctuations (see the latter is small enough 
to be covered in the data points for average). Interestingly, the average velocity of the 
numerical data for fS (diamonds) are also on the curve in the other side, even though 
there are considerable fluctuations as indicated by the bars. This suggests fS can also be 
understood with HVS. Below, we demonstrate that the dynamic property of HVS can explain 
this observations. 

We performed the linear stability analysis 18] on HVS and find this is linearly stable, 
including marginal stability, regardless of density and model parameters (see Appendix B). 
When v/rr] > 1, there are at most two marginally stable modes out of the total 2 N stable 
eigenmodes, where N is the number of vehicles. Otherwise with v/rr] < 1, a half of the 
total modes are marginally stable. Thus when HVS is realized with N 1 with v/rr] > 1, 
it readily shows almost uniform velocity over the whole system while, with v/tt] < 1 , 
it may exhibit the fluctuations attributed to the macroscopic number of marginal modes. 
This dynamic property is consistent with the numerical observation on hS and fS shown in 
Fig. EKa) 19]. 

ft is worthy of noting that the stability boundary of v/rr] = 1 is identical to that of the 
numerical phase boundary shown in Figs. [2](b) and[3](a). When v is replaced with p using 
Eq. ()6l) , (p ^ 1 — L)/r 2 r] = 1 is immediate. This is the phase boundary (solid curve) shown 
in Fig. |2](b) , for r = 1. The macroscopic number of marginal modes in fS can explain the 
observation of J in the fS-region (see Fig. 2(b)). Since the findings so far hold for any r 













(1/p - L)/t 2 t| 


FIG. 3. (Color online) Density-velocity relation and deceleration, (a) Scaled density-velocity 
relation for mu = 5 m/s 2 and various r = 1, 1/2, 1/4, 1/8, 1/16 sec. The solid line is the analytic 
curve by Eq. ([6]). Each point is the 10000-sec. average of vehicles’ velocities in steady state for a 
given density (all the points in Figs. 3 and 4 are obtained in the same way). The bar shows the 
fluctuation of the velocity, (b) Average, rms (root-mean-square) fluctuation, and maximum of the 
decelerations in the data giving (a) (each deceleration is measured in the ratio to the deceleration 
capacity). 

and statistics of D n s, we conclude that hS and fS are the dynamic phases of our model. We 
finally emphasize the phase boundary condition of v/rr] = 1 is same to the condition where 
at least one vehicle follows u d . The flow established in the presence of such a vehicle is the 
very hS that is much more stable than fS. This indicates that u d unrecognized in the earlier 
models plays a significant role in stabilizing the whole system. 


V. COMFORT AND FLUX 


In the following, we examine our model generates a practically appealing traffic flow. If 
traffic flow is safe, one of the next concerns is the driving comfort, which is required for 
autonomous driving systems 


2Q 


2l| . For this, we measure the decelerations each vehicles 


experience in the simulation for Fig. [3^a). Each deceleration is measured in the ratio to 
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FIG. 4. (Color online) Effect of response time, (a) Average speed as a function of the density p 
and the response time r when mo = 5 m/s 2 . The plateau in the upper left corner is speed limit 
Umax. The solid curve separating the data point type is the phase boundary between hS and fS, 
and the dashed curve below is its projection onto the p-r plane, (b) The projection of (a) onto the 
density-flux plane for t = 1,1/2,1/4,1/8, and 1/64 sec. The upper-bounding curve is obtained in 
the r —} 0 limit [see Eq. ©]. All curves in (a) and (b) are analytic results. 


the deceleration capacity. The results are shown in Fig. (3][b). We obtain three statistical 
observables; average, fluctuation (root-mean-square), and maximum of the ratios. The 
average is approximately 0.1 and 0.001 ~ 0.01 in hS and fS flows, respectively, and the 
root-mean-square shows the similar values. The maximum is around 0.01 ~ 0.1 and 0.2 ~ 
0.7, respectively. We remark the deceleration strength is a characteristics of flow phase as 
observed, and thus driving comfort can be considerably improved by promoting hS with 
smaller r (see Fig. [4]). 

The other practical interest is probably the flow efficiency, which can be represented by 
vehicular flux. For a homogeneous-velocity solution v, the flux is pv by the hydrodynamic 
relation jyj]. Since 2 f(z) > zf\z ) as a general property of f{z) (Appendix A), the steady- 
state velocity v increases as r decreases. The velocity-increase for smaller r is drawn in 
Fig. [U(a) along each constant-/? curves. This gives the flux-increase, as shown in Fig. [4](b). 
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Since f(z) ~ z 2 for large z, the flux converges to 


cVp(i/£-p) 


( 7 ) 


in the r —> 0 limit [see the upper-bounding solid curve in Fig. [4(b)], where C is a constant by 
the statistics of D n s [see Eq. (1A26h for the detail]. We consider this result is also appealing 
because i) the flux for r = 1 sec. (the typical response time of drivers 22]) is comparable 
to empirical maximum value around 2500 veh/h 23-4251]. ii) the flux enhancement is more 


sensitive for larger r, iii) the flux becomes considerable for r of O.l-sec.-order, and iv) all 
these are achieved in the manageable condition of Eq. ([5]) guaranteeing safety. It is worthy 
of noting that the research held of autonomous driving system has already been treating the 


processing time down to 0.1 sec. 20 


211 • 


VI. FINAL REMARK 

We finally remark that an extension of our model in order to cover the other features of 
traffic how (lane change, on-ramp how, traffic signal, and so forth) is straightforward. This 
is because the problem is still a compromise between safety and mobility in consideration of 
the positions, velocities, and deceleration capacities of the related objects. Also, our model 
provides a new perspective to the study of traffic-how safety through the interpretation of 
the non-manageable events and its statistics. Besides, we expect an autonomous driving 
system based on our model can be possible in the solid safety criterion and in the appealing 
traffic-how quality of comfort and hux. 
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Appendix A: Homogeneous Solution 


Let a velocity v be the homogeneous solution of our model. Substituting it for all the 
velocities in Eq. ([2]), we obtain the optimal spacing as 


TV 


S n = 


D 


71 + 1 


D„ 


1 _3_ 

^n+l D n 


for 

for 


l 


D n +i 

1 

Dn + l 


—< I 
D„. — v’ 


— > - 

D n V 


(Al) 


Then the average spacing is S — Y2n =i A^/A’, where A is the total number of the vehicles. 
The (global) vehicular density is simply given by 

1 


P = 


S + L ’ 


and the average flux is 


q = pv. 


(A2) 


(A3) 


Equation (1A11) can be written as 

, . s s (v, y) — tv — \yv 2 for y < t/v, 

s(v , 2 /) = % 2 (A4) 

y) = ij for y > t/v , 

where y stands for A d n = l/D n+ i~ 1/D n . When D n is randomly assigned out of the interval 
{me — wd, m D + wd), the range of A d n is (—1 /rj, 1 / 77 ) where 


1 _ 1 1 
r] m D — w D m D + w D ' 


(AS) 


Note that u d (and accordingly s d ) does not appear when v < ti 7 . 

If there are sufficiently many vehicles and their D n ' s are uncorrelated, the average spacing 
can be obtained by the integral 

f 1 A f T / v ( v 2 y\ f 1/v t 2 

S= s(v, y) P Ad (y) dy= ( tv -— ) P Ad (y) dy+ —P Ad (y) dy, (A6) 

J — l/r) J —l/rj \ " / Jr/v ^V 

where P Ad (y ) is the probability density function of A d n = 1 /D n+ i — 1/D n . Thus the average 
spacing of homogeneous solution with v is completely determined by the distribution of A d n , 
which can be obtained directly from the distribution of D n . Let us denote the probability 
density function of D n as Po{y)- Then P d (y) of 1 / D n is given by P d (y) = ( l/y 2 )P D (l/y ) 
and P Ad (y ) is given by the convolution P Ad (y) = P d {x) P d (x — y) dx. Note that PA d (y) 
is an even function, P Ad (-y) = P Ad (y). 
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Let us introduce the dimensionless speed z = v/ry and the rescaled probability function 
p(y) for rj(l/D n+ i — 1 /D n ). Then Eq. (1A6|) can be written as 

£; = L' j l/z T v v{ ' y)dy - (A7) 

Since p(y) has the normalization j\p(y) dy — 1 and the symmetry property p(—y) = p(y), 
we have f^ z p(y) dy = 1 — J^_p(y)dy and f// z yp(y) dy = — p(y) dy. Therefore we 
obtain 


F = z[l - M*)] + jU*) + \l- 1W = F(Z), 


where 


h{z)= / y k p(y) dy. 

Jl/z 

This provides the scaling relation between S (or p) and v, 


x = 


S 

T 2 rj 


P 


-l 


T 2 T] 


= F 


TT] 


= F(z). 


(AS) 

(A9) 

(A 10 ) 


Now we find the general properties of F(z). Since p(y) = 0 for \y\ > 1, we obtain 
Ik(z) = 0 for z < 1, leading to F(z) — z for z < 1, or equivalently, q — (1 — pL)/r for 
P > 1 /(t -2 77 + L). Therefore we arrive at 


F(z) = 


for z < 1 , 
f(z) for z > 1 , 


(All) 


where the function f(z) is determined by the distribution p(y). Since h{z = 1) = 0, we 
obtain lim z4l + f(z) = 1 = F(l). Thus F(z) is a continuous function. 

We have f(z) = z/2+ f Q 1/z z + j\ /z {z 2 y /2 + 1/2 y) p(y) dy from I 0 (z) = 1/2-f Q 1/z p(y) dy. 
Since Jq p(y) = 1/2 and the integrand function has the minimum z and the maximum 
(. z 2 + l)/2 over [0, 1], the integral is bounded between z/2 and ( z 2 + l)/4. Therefore we 
obtain 

z<F(z)< \(z+l) 2 . (A12) 

The lower bound F(z) > z leads to v < (p -1 — L)/t, which gives the upper bound of the 
average flux q < (1 — pL)/r. 

For z < 1, we have F\z) = 1. Using 


I'k{ z ) = ^J ykp (^) dy = 


Z k+ 2P U 


(A13) 
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we obtain f'[z ) = 1 - I 0 (z) + z/i(z) = 1/2 + / 0 1/2 1 + J^ z zy P(y) dy. Since the integrand 
has the minimum 1 and the maximum 0 for z > 1 , the integral is bounded between 1/2 and 
z/ 2. Thus we have 


1 < f(z) < 


z + 1 


(A 14) 


Since F'(z ) > 0 for all z > 0, there exists the inverse function G = F 1 such that G(F(z)) = 
G(x) = z. Therefore we obtain the inverse relation of (1A10|) as 


z = — — G 

TT] 


P 


-1 


T 2 7] 


L\ I x 

= G{x) = 


for x < 1, 
g(x) for x > 1. 


(A15) 


From (1A12I) and (IA14F we obtain 2y/x — 1 < G(x) < x. Since G'(x) = 1 /F'(z), Eq. (1A14h 
leads to \fx < g'(x) < 1. From dv/dp = —G'(x)/rp 2 } we obtain —1/rp 2 < dv/dp < 0. 
From (IA1 jj) and (1A131) . we obtain 


F"(z) = h(z) - I' 0 (z) + zl[(z) = h(z) > 0, 


(A16) 


leading to G"(x) = — F"(z)/(F'(z)) 3 < 0. This implies the average flux q is a non-concave 
function of p because d 2 q/dp 2 = G"(x)/r 3 r]p 3 < 0. 

From (1A15[1 . dv/dr at a fixed p is given by r)(G(x)— 2xG'{x)) = —rj(2F(z)—zF'{z))/F'(z). 
For z < 1, we have 2 F{z)—zF'{z) = z > 0. From (1A6|) and (1A14I) . we obtain 2 f{z)—zf'{z) = 


(Air) 


z (1 — Io(z)) + I-i(z) = z + z fy z (1 /zy — 1 )p(y) dy > z. Thus we have 


2 F(z) - zF'(z) 
n <_Ed_< z 

U < F'(z) S ’ 


leading to — (p _1 — L)/r 2 < dv/dr < 0. 

There is another lower and upper bound of F(z). Using Jo(oo) = 1/2, we have f(z) = 
z 2 1\ (oo )/2 + z/2 + Jq / ' Z (z — z 2 y/ 2 ) + f^ z (l/2y) p(y) dy. Since the integrand has the max¬ 
imum x and the minimum 1/2, the integral is bounded between 1/4 and z/2. Therefore we 
obtain 

f 1 r 2 

(A 18) 


- 1 -< f(z )-ii(oo) < 2 ). 

2 4 w 2 n ~ 


This leads to the asymptotic relation 


F(z) ~ A_/i(oo) for 1, 


(A19) 


which implies q ~ y/ (2r]/Ii(oo))p(l — pL ) for p 1 -L> r 2 p. 
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If D n has uniform distribution on {mo ~ u>d, m d + u>d) = (a, b), then Pd(u) = 1/2 Wd, 
which leads to Pd(y) = 1 / 2 wr>y 2 . Thus we obtain 


PAd(y) = 


R(\y\) for \y\ < 1/r 7 , 


0 


where 


R(y) = 


for \y\ > 1/r], 
b ~ a ~ + rfe ~ l ln [(1 ~ ay)(l + by)} 


4 w 2 D y 2 


(A 20 ) 


(A 21 ) 


After some algebra, we obtain 

1 f 2 


/(*) = z 


7]Z 


y 2 z 3 


iwl 1 3 (1 - z) ( “ te + - 3 


111 | 1-1 + ln (1 H- 

yz J \ yz 


, , b b \ fa a 

a In [-- b In - 4- - 

a yz J \b yz J j or/ |_ \a y 


a 3 ln f— — — ) — b 3 In I 2-1 + - 


az a 


b y 


(A 22 ) 


In the limit of r —> 0, we have s(v,y ) = —^yv 2 for y < 0 and 0 for y > 0. Then the 


average spacing is given by 


where 


St —>0 ’ 




2y 


fo yp(y)dy 


(A23) 

(A24) 


Then the average speed and flux are given by 


V (p _ L )’ Qt ^° = C ^ p(yl ~ pL ^' 


'Vt—> o c j 

For the uniform distribution of D n over {mo — Wd,^d + Wn), we obtain 

C 2 = 


(A25) 


2 w 2 d 


m D tanh 1 ( — 
u V m D 


(A26) 


w D 


Appendix B: Linear Stability Analysis 


We investigate the linear stability of the homogeneous solution with respect to small 
perturbations. Assuming that the velocity and spacing are very close to those of the homo¬ 
geneous solution, we can write down 

v n (t) = v + u n (t), s n (t) = S n + a n (t), (Bl) 
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where the optimal spacing S n is given by (lAlj) . By linearizing (3), we obtain 

a n(t) — 2 Un (0 D n+ ^n+lU) 


v s n(t + r) = v + 


V n( t + T ) =V + U n+1 (t) + 


T I V 

2 ^ D n 

D r, - D 


n+1 


D n + D 


n+1 


2cr n (t) 


^n(^) “ 1 “ ^n+]+) 


(B2) 


up to the first order of u n and a n . Meanwhile, from the integration scheme (1), we obtain 

T 

°n{t + t) = a n (t ) + - [u n+1 (t) + u n+1 (t + r) - Unit) - Unit + t)] . (B3) 

The periodic boundary condition Xn+i = X\ is applied for any quantity x. 

We introduce the new variable 


I /,n Cnif) U n (t) U n +i(t) 

) = — 2 2 ' 

From (IB3D . we obtain 

+ r) = 1p n (t) ~ Unit) + U n+1 (t). 

Then we can rewrite (1B2I) as 

s /, \ ++) u n{t) OL n +1 /,\ 

v n (t + r) - v = -1- u n+ i{t), 


0L r . 


Oi r , 


Oir 


1>J(t + t) - V = Pnii n (t) - P n U n (t) + (1 + p n )u„ +1 (t), 


(B4) 


(B5) 


(B6) 


where 


Q+ — 


+ = 2 


Dn — D 


n+1 


(B7) 


2 rD n D n + D n +1 

Note that a n > 1/2 and 0 < l/a n < 2. The selection criterion 1 /D n+1 — 1 /D n > r/v for + 
is equivalent to a n +i — a n > 1. Now we combine (IB 01) into 


where 


and 


Un(t + t) = 7 n^n(t) ~ 7 + #n U„ + i (t), 

In = 7n) + = #n for «n+l ~ Oi n < 1, 

7n = 7n> ^ for «n+l “ > 1, 


s = + _ = 

" «n’ 71 OL n 

7n = Pn, °n = 1 + + 


(B8) 


(B9) 


(BIO) 
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Here we have assumed that the selection of v s n or v ^ is not changed by u n . Note that 7 „ and 
9 n are continuous functions of a n and a n+ 

By introducing the vector notation f(£) = (?/q(f), ■ ■ ■ , ^jv(f)), u(i) = • • • , Wjv(t)), 

and x(t) = (f (f), u(£)), we can combine (IB5D and (IB8D into a Jacobian matrix equation 


x(£ + r) 


( f(* + r)\ 
\u(t + t)J 


fl N B\ /f(t)\ 
V C T ) \n(t)) 


J x(t), 


(Bll) 


where Ijy is the N x N identity matrix and the lower-left submatrix C is a diagonal matrix 
with G nn = 7 „. The upper-right sub matrix B and the lower-right submatrix T are given by 



f-1 

1 

0 •• 

• 0 N 


^-71 

9i 

0 •• 

• 0 ^ 



0 

-1 

1 • • 

• 0 


0 

-72 

e 2 •• 

• 0 


B = 

0 

0 

-1 •• 

• 0 

T = 

0 

0 

-73 •• 

• 0 

(B12) 



0 

0 •• 

• -V 


V 9n 

0 

0 •• 

• -In) 



The long-time behavior of the perturbation amplitude is determined by the largest mag¬ 
nitude among the eigenvalues. The eigenvalues are given by the zeros of the characteristic 
polynomial 


h{ A) = det (J — AI 2 tv) = det 


A B 
C D 


= 0, 


(B13) 


where A = (1 — A)Ijv and D = T — AIjv- Using the properties of the determinant, we obtain 


N 


N 


h(\) =a k n ii n (\) - n u n (\)= o, 


(B14) 


n= 1 


n= 1 


where 


/^n( / ^) — ^ 1 “h Tm ^n(^) — @n 1) “t“ 0V? 


(B15) 


For A = 1, we have /i n ( 1) = cu n (l) = 7 n , leading to h( 1) = 0. Therefore at least one 
eigenvalue is exactly 1. From (1B14D , we have ti(l) = (Iln=i 7 n) Z)«=i ( 1 + (1 - Qn)/ln)- 
Since the summand is zero for v ^ and positive for u®, we obtain h'( 1) > 0. Thus the 
eigenvalue 1 is unique without degeneracy. On the other hand, for A = —1, we have 
co n (— 1 )/ fi n (— 1 ) = ( 2 a n+ i — l)/( 2 a n — 1 ), irrespective of whether w® el is u® or leading to 
/ fi n (-l)) = 1. Thus we have h(- 1) = (1 - (-1) A ) E[f=i( 2 _ 7 n), which is zero 
for even N and negative for odd N. Thus —1 is an eigenvalue if and only if N is even. For 
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even N, we obtain h'(-l) = (-l) N {Un=i ( 2 ~7n))(~N + En=i ( 1 /( 2 ~ln/Q n ) ~ V ( 2 ~ 7 n)))- 
The sum is zero when u® el = v s n for all n. Since 1/(2 — 7 ^/^) — 1/(2 — 7 ^) < 1 /(2 — 7® / 9 s n ) — 
1/(2 — 7 ®) for a n+ \ — a n > 1, the sum is not positive. Thus we obtain h!{— 1) 7 ^ 0 and —1 
is also a unique eigenvalue. 

Moreover, ±1 are the absolute bounds of real eigenvalues. For A > 1, we have /i n (A) > 0 
and u; n (A) > 0. Similarly, we have /z n (A) < 0 and u; n (A) < 0 for A < —1 since 7 „ — 2 9 n < 0. 
For A to be a solution of (IBldjh it is necessary to satisfy | j/[^ = 1 u; n (A)//i n (A)| = |A|^. 
Meanwhile, it can be easily shown that 

N 




n 

n= 1 


<4(A) 


<(A) 


= 1 , 


(B16) 


due to the transitivity a;® (A) = 0 s n n s n+ X (A) and = 1- For |A| > 1, we obtain 

^n(A) ^ w®(A) 


/4(A) 


< 


/4(A) 


for ot n -)-i 7 s 1 , 


(B17) 


by using 


KWrf(A)|-K(A)<(A)| = 


(^n+1 ^n)(^n+1 




1 ) 


(|A| 2 -1). (B18) 


®n(®n+l T 1) 

Therefore we obtain n^Li |ca n (A)//i n (A)| < 1 < |A|^ for |A| > 1. Thus there is no real 
eigenvalues in the range A>lorA<—1. 

The equality (IB 161) and the inequality (IB 171) holds even for the complex A = re ^ with 
the magnitude r > 1 and angle 0 (0 < 0 < 2tt). After some algebra, we have 


k(re‘*)rf(r e «)| 


\ui(re‘*) A(re‘t)\ 

(®n+1 C^n)[(c^n+1 ®n) 1] 


a n( a n +1 + Otn ~~ 1) 

= M n (r, COS 0), 


4 (l — cos 2 0) + (do + ai cos 0) (r 2 — l) + 02 (r 2 — l) ‘ 


(B19) 


where a 0 , cq, and a 2 are constants to be determined by a n and a n+ i. Due to (IB 171) . we have 
M n (r, ±1) > 0 for r > 1 , which corresponds to A = ±r (0 = 0 or n). Since M n (r, cos0) is 
a quadratic function of cos0 with the negative coefficient on the quadratic term. Therefore 
M n {r > 1, cos 0) is positive over the entire range of 0. This leads to 

a; d( re «0) cu®(re^ 


/4( rei 


< 


/4/re^) 


for OL n .y\ Ol n k 1, 


(B20) 
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for any angle 0 if r > 1. Thus we can conclude that |^n( r 'C 9i )//i n (re* 9;, )| < 1 < r N for 

r > 1. Therefore there is no complex eigenvalues in the region |A| > 1. 

If v d does not appear at all (p > l/(T 2 r] + L)), we have = l/a n and 9 n = a n+ i/a n 
for all n. Then we have p n (X) = X — 1 + l/a n and cu n (A) = (a n+ i/a n )(A — 1) + l/a n = 
a n+ ip n+ i(X)/a n , leading to nlili^A) = n»=i/ i 2 3 4 n(A). Thus we obtain 


N 


n= 1 


(A"-l)n A-l + -Uo, 


(B21) 


which leads to A = e 2mn ^ N or A = 1 — l/a n (n — 1, • • • , N). Since |1 — l/a n \ < 1, we obtain 
| A | m ax = 1- Thus a u s -only homogeneous flow is marginally stable. Note that N eigenvalues 
among the total 2N eigenvalues have |A| = 1. 

On the other hand, if u d is selected for at least one vehicle, the situation changes drasti¬ 
cally. For A = e*^, we obtain 


k(^)rf(^)r - kkKk)f = 


2 (c^n+1 C^n) [(otn+1 1] 


a n( a n +1 + Oi n ~ 1 ) 


4(1- 


cos 


(B22) 

by substituting r = 1 in Eq. (1B19I) . Since — 1 < cos0 < 1 for any non-real A along the unit 
circle on the complex plane, we have 


(n d (e^) 


/4(e^) 


< 


<(e^) 


/4(e^) 


for cx n 1 , 


(B23) 


unless 0 is an integer multiple of 7 r. Thus we have n»=i \ u n{^) / 'p n {^)\ < 1 for all non-real 
A = e ^ if at least one v d is selected. Therefore, all the eigenvalues except for 1 (and — 1 
for even N) have the magnitude less than 1. Consequently the w d -mixed flow is much more 
stable than the u s -only flow. 
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